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Abstract: We present observational results from optical long-slit spectroscopy of parsec-scale jets of DG 
Tau. From HH 158 and HH 702, the optical emission lines of Ha, [O i] AA6300, 6363, [N n] AA6548, 
6584, and [S n] AA6716, 6731 are obtained. The kinematics and physical properties (i.e., electron density, 
electron temperature, ionization fraction, and mass-loss rate) are investigated along the blueshifted jet 
up to 650" distance from the source. For HH 158, the radial velocity ranges from —50 to —250 km s -1 . 
The proper motion of the knots is 07196 — 07272 yr _1 . The electron density is ^10 4 cm -3 close to the 
star, and decreases to ~10 2 cm -3 at 14" away from the star. Ionization fraction indicates that the gas 
is almost neutral in the vicinity of the source. It increases up to over 0.4 along the distance. HH 702 is 
located at 650" from the source. It shows ~ —80 km s _1 in the radial velocity. Its line ratios are similar 
to those at knot C of HH 158. The mass-loss rate is estimated to be about ~ 10 - ' M 0 yr _1 , which is 
similar to values obtained from previous studies. 

Key words: stars: formation — stars: outflows Herbig-Haro obejcts — stars: individual: DG Tau 
techniques: spectroscopy 


1. Introduction 


In the formation of stars and planetary systems, 
studying outflows a nd mass-accretion is essential 
(lHartigan et al.l H995I ). Through the analysis of the 
kinematics and physical conditions of outflows from 
young stellar objects (YSOs), we are able to understand 
the jet launching mechanism and the interaction with 
ambient material. Parsec-scale jets are useful to study 
star formation in somewhat long time scale and large 
spatial distribut ion. There are a number of parsec-scale 
jets known (e.g., Ravlll987l : lBallv fc Devinelll994IOgural 


fl995t iBallv et all 199. r jj ~ One of the extreme cases is the 

HH 222 system, suggested as a giant Herbig-Haro flow 
s pread over 5.3 pa r sec in length rtReinurth et al.ll2013l) . 

iMundt fe Friedl (119831) discovered a jet-like outflow 
from DG Tau. DG Tau is an early evolutional stage 


CTTS (Class II) with mass of 0.67 M@ (Hartigaii£t_al 


1995 ) and age of about 3 x 10 5 yr dBeckwith et al 


1990). IMcGroartv fc Ravi ( 20041) suggested that it is 


the driving source of a parsec-scale jet including HH 
158, HH 702, and HH 830. The later two HH objects 

are located at more than 9' away from DG T au. _ In 

their follow up observations, IMcGroartv et all 12 0071 ) 
reported that HH 702 showed receding motion with re¬ 
spect to DG Tau. They also concluded that HH 830 is 
not ejected from DG Tau. lEisloffel fc Mundtl ( 19981) 
showed that the HH 158 outflow extends to ~ 11" 
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toward the southwest direction and estimated the di¬ 
rection of the proper motion of knots from DG Tau 
as ~226 ± 7°. They also calculated the jet inclina¬ 
tion angle as ~ 32° with resp ect to the line o f sight . 
In addition to an opti cal jet, iRodrfguez et al.l (2012) 
and iLvnch et al.l ( 20131) detected a radio knot at ~ 7" 
from the central source. High spatial resolution obser¬ 
vations with a space telescope and ground based large 
telescope with an adaptive optics system showed that 
the blueshifted H H 158 jet can be traced to within ~ 
071 from the star (iBacciotti et al.ll2000l 20021 : iPvo et all 


120031 iMaurri et ah l2014lh Lavallev et alJ ( 1997 ) re¬ 
ported^ that the redshifted jet with the length of ~ 1". 
IPvo et akl (120031 ) showed that there is a 077 gap be¬ 
tween the redshifted jet and the star in the [Fe ii] 1.644 
/im emission, which indicates the presence of an op¬ 
tically thick disk. The two velocity components are 
shown in t he optical and near infrared forbidden emis¬ 
sion l ines (Hartm ann fe Ravmondl Il 989: Solf & Bohm 


1993; 


2003, 


ncs (Hartmann Raymond 1989; son & Eomr 

Hamann 1994 Hartigan et akl 19951 Pvo et al 
120061 ) . 


In this paper, we report the results of optical long- 
slit spectroscopy on HH 158 and HH 702 outflows. We 
investigate the electron density, electron temperature, 
and ionization fraction from close vicinity of the central 
source to the knot at 14" of HH 158 jet and HH 702 at 
650" from DG Tau. We note that this is the first time 
to obtain optical spectroscopy of the HH 702 outflow. 
We also summarize the proper motions of the knots in 
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Table 1 

Log of Observations. 


Object 

Date (UT) 

P.A. O 

Exposure 
time (s) 

Slit 

width (") 

Spectral 

resolution (A/A A) 

Spectral 
coverage (A) 

HH 158 

2012 Nov. 15 

223 

1800 

2.9 

2070 

5900-7050 

HH 702 

2012 Nov. 15 

193 

2800 

2.9 

2070 

5900-7050 



Figure 1. Slit positions on the object field. The center of 
the 2'.'9(W) x 3l6(L) slit is set on HH 158 and HH 702. 
The position angles are 223° and 193° for HH 158 and 
HH 702, respectively. The background image is from the 
Second Palomar Observatory Sky Survey (POSS II). The 
HST image of HH 158, 1 59 and Ha image of HH 702 from 
IMcGroartv fe Ravi (120041 1 are inserted as insets. Knot A — 
E of HH 702 are marked. 


the jet with the positions during the last 30 years. In 
Section 2, we describe the observation and data reduc¬ 
tion. In Section 3, we show the line intensity and ratios 
in the position-velocity map. Sections 4 and 5 are the 
discussion and summary, respectively. 

2. Observation and Data Reduction 

The observation was conducted on November 15, 2012 
(UT) at BOAO (Boliyunsan Optical Astronomy Ob¬ 
servatory) with an optical long-slit spectrograph in¬ 
stalled on the Cassegrain focus of a 1.8-m telescope. 
The wavelength coverage was 5900 — 7050 A with 1200 
gmm -1 grisrn, which was chosen to achieve the emis¬ 
sion lines of Ha, [S ii] AA6716,6731, [N n] AA6548,6584, 
[O i] AA6300,6363. The slit length and width were 3.6 
arc-minute and 2.9 arc-second, respectively. The re¬ 
sultant spectral resolution was ~ 2070 (150 km s -1 ) 
with the dispersion of 0.41 Apixel -1 . The CCD cam¬ 
era size is 4k x 4k and the pixel scale along the slit 
is 0'.'45 pixel -1 . Exposures on comparison lines were 
taken with the FeNeArHe lamp before and after the 


DG Tau ID Spectrum 



Figure 2. Integrated ID spectrum of the central region of 
DG Tau. 


exposure of each object, for wavelength calibration. 

Table [I] shows the observation log of HH 158 and HH 
702. Figure [1] shows the slit positions. The position 
angle (PA) of the HH 158 re gion including the centra l 
star is determined to be 223° (IMcGroartv fe Ravll2004[) 
after verification through an HST image. In the case 
of HH 702 located at ~ 11' from the star, a PA of 
193° was selected to obtain the emissions as much as 
possible. The exposure times were 1800 s and 2800 s 
for HH 158 and HH 702, respectively. HR 1544 (Al 
type) was observed for telluric line correction and flux 
calibration. 

We reduced the data with standard IRAF (Image Re¬ 
duction and Analysis Facility) packages. Bias was sub¬ 
tracted from every frame. The flat-fielding was carried 
out using a normalized flat made by the APFLATTEN 
task. Bad pixels and cosmic-ray events were corrected. 
The spectra were extracted using the APALL task. The 
wavelength calibration and distortion correction were 
accomplished by the IDENTIFY, REIDENTIFY, FIT- 
COORDS, TRANSFORM tasks. There were 1 and 2 
pixel shifts between the comparison frames before and 
after the exposure for HH 158 and HH 702, respectively. 
Since the instrument is mounted on the Cassegrain fo¬ 
cus, the instrumental displacement due to changes of 
the telescope elevation angle cause those shifts. We 
used the average of two lamp images to achieve a more 
accurate spectral calibration. There were also wave¬ 
length shifts between object frames. We corrected these 
shifts with the OH airglow lines in the object frames. 
After that, the sky emission lines were subtracted using 
the BACKGROUND task. The wavelength sensitivity 
correction and flux calibration were done with the stan¬ 
dard star, HR 1544 (Al, V=4.35, T eff =9150 K). 
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Figure 3. Continuum-subtracted PVDs of the emission lines of HH 158. Contours are drawn from 0.001 (3 cr) to the peak 
levels of 5.2 (Ha), 0.98 ([O i] A6300), 0.33 ([O i] A6363), 0.17 ([N n] A6548 and A6584), 0.10 ([S n] A6716) and 0.22 
([S n] A6731) x 10~ 15 W m -2 A -1 with an equal interval in a logarithmic scale. The position angle (223°) of the slit 
is marked in the left panel. Northeast is up and sout hwest is down in the diagram . The original spectrum has Gaussian 
Smoothed with a a value of 1.5 x 1.5. Knots Bl, B2 llLavallev et al .1 Fl 997l ) . and C (lEisloffel fc Mundtll 19981 1 are identified 
at Y ~ —6” and —14”. The faint emissions at Y = —54” in [O i] A6300, [S n] A6716 and [S n] A6731 are parts of the HH 
159 outflow ejected from DG Tau B. These are labeled with the ‘x’ symbols in the figure. 



Figure 4. PVDs of the emission lines of HH 702. Contours are drawn from 0.05 to the peak levels of 1.5 (Ha), 0.45 
([O i] A6300), 0.12 ([O i] A6363), 0.11 ([N n] A6548), 0.38 ([N n] A6584), 0.60 ([S n] A6716), and 0.58 ([S n] A6731) in 10" 17 
W m -2 A -1 with an equal interval in a logarithmic scale. The position angle ( 193°) of the slit is marked in the left panel. 
Northeast is up and southwest is down in the diagram. The knots A and E of lMcGroartv fe Ravi (|2004l l are identified. 

























































116 


Oh et al. 



Figure 5. (a) Variation of the line profile, ( b ) peak intensity, 
(c) peak velocity, and ( d ) velocity width (FWHM) of the 
each emission line of HH 158 with the distance from DG 
Tau. 

We removed a deep Ha absorption line in the spec¬ 
trum of the standard star using the SPLOT task before 
applying the wavelength sensitivity correction. 

3. Results 
3.1. Kinematics 

Figure [5] shows the ID spectrum extracted from the 
central source region within ±5"4 from DG Tau. In 
the wavelength range of 6200 — 6800A, the Ha A6563, 
[O i] AA6300, 6363, [N ii] AA6548, 6584, and 
[S n] AA6716, 6731 emissions were detected. The 
[Fe II] emissions at A6433, A6456, A6516 A and the 
atomic Helium line at A6678 A are also identified, the 
Ha emission has the highest intensity peak level of 6 
x 10 -14 W m -2 A -1 , which is ~ 30 times higher than 
that of the stellar continuum. The wing of the strong 
Ha line affects the line profiles of two [N ii] lines. 

Figures [3] and 0] show the position-velocity diagrams 
(PVDs) of HH 158 and HH 702. In Figured the stellar 
continuum was subtracted. The position at Y = 0" 
indicates DG Tau and the peak of knot E in Figures 
[3] and 21 respectively. HH 702 is located at ~ 11' (= 
0.45 pc at 140 pc distance) southwest from DG Tau, 
as shown in Figure [TJ The Gaussian smoothing was 
applied in both figures with sigma of 1.5 x 1.5 pixels. 

In Figure p w e identi fy the knots Bl, B2 
dLavallev et allfl997h and C feisloffel fe Mundtlll998l) 
at the distance of 6'.'74, 8'.'5 and ~ 14" from the star. 
The two velocity components appear in all emission 
lines in Figure [3j The low velocity component (LVC) 
is located within ±5". The peak velocity of the LVCs 
varies from ~ —80 to —20 km s -1 . The high veloc¬ 
ity component (HVC) is relevant at |Y| > 6". Fig¬ 
ure [5]c shows that the HVC is getting slower with in¬ 
creasing distance, e.g., —270 km s” 1 at 6"3 (knot Bl) 
and —150 km s _1 at 18"45. The real jet velocity is 
1.18 times higher than the radial velocity considering 



V l sr (km s' 1 ) 


Figure 6. Line profiles at knots A (upper) and E (lower) of 
HH 702. 

the in clination angle of 32° wi th respect to the line of 
sight (iEisloffel fe Mundtl Il998lh so the radial velocity 
of —270 km s _1 corresponds to the jet velocity of —320 
km s . Additionally, it is noticeable that the emission 
lines are showing different spatial extents. In Figure [3l 
all emissions except [N ii] A6548 show extensions to Y 
~ —20" at 3 a contour levels. Knot Bl is relevant in all 
emissions except [O i] A6363. Knot C is remarkable in 
[S ii] A6731 and appears as plateaus in Ha, [O i] A6300, 
[N ii] A6584, and [S ii] A6716. Its intensity in [N ii] 
A6584 is almost twice than that of [O i] A6300. 

The faint emissions at a distance of 54" from DG Tau, 
detected at [O i] A6300, [S ii] A6716 and [S ii] A6731 
may not come from HH 158 but from a part of HH 159 
(see Figure [T]), which is ejected from DG Tau B. These 
features are labeled with ‘x’ symbols. 

In Figure IH we ident ified the knots A and E of 
iMcGroartv fcRavl (|2004f ). The knot A of HH 702 is 
bright only in Ha emission, wh i le it is s een in the 
[S_ii] images of lMcGroartv fe Ravi (l2004h and 1 Sun et al.l 
(l2003f) . Knot E is detected in all emission lines. 

Figure [5] shows the variations of the line profile, peak 
intensity, peak velocity, and velocity width (FWHM) 
of each emission line with the distance from the star 
for HH 158. The data are sampled at 6"30, 10"35, 
14"40, and 18"45 along the blueshifted jet. The plot 
starts from 6"30 where the knot Bl is located. The 
[O i] A6300 emission at 18"45 could not be measured 
due to faintness. In Figures [5}5 and [5]s, the peak in¬ 
tensities and velocities of all the emission lines decrease 
with distance from the source. The uncertainty in the 
peak velocity measurement of [N ii] and [S ii] at 18"45 
is ±170 km s _1 . It is large, due to low a signal-to-noise 
ratio (S/N < 2). Even though there are small increases 
in velocity (Figure [5]c) , those are not significant due to 
a large uncertainty. The decreasing rate in the velocity 
of Ha is \Av peak \ ~7 km/". 

In Figure [5ji, the range of velocity widths (FWHM) 
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Figure 7. Ratios of forbidden emission lines: (top) HH 158 and ( bottom ) knot E of HH 702. (Left) [S n] A6731/ A6716, 
(center) [N n] A6548+A6584 / [O i] A6300+A6363, and (right) [O i] A6300+A6363 / [S n] A6716+A6731. 


varies from 150 to 350 km s — , which is marginally re¬ 
solved considering the velocity resolution of 150 km s _1 
in our observation. The low velocity component shows a 
wider width (~340 km s _1 in Ha) than those estimated 
for high velocity components (~160 — 230 km s _1 ), 
comp atible with results obtained i n pre vio us studies 
(e.g., iHirth et al.1 Il997t iPvo et all l2003h . iPvo et all 
( 200211 interpreted the wider velocity width of LVC as 
follows: if we assume that the outflow has the shape 
of a diverging wind, the velocity width represents the 
difference in radial velocities between the nearest and 
farthest streamlines from us. Thus, the wider velocity 
width means a wider opening angle. 

Figure [S] shows the variations of the line profile along 
the distance of HH 702. The peak velocities of two 
knots are ~ —80 km s -1 in the Ha emission. 

3.2. Line Ratios 

The diagnostics by the ratios of forbidden emission 
lines is a very effective tool to study .jet properties 
(lHartigan et al.|[T994 iBacciotti et al.lll995lh In Figure 
0 we show the flux ratios of the forbidden emisson lines 
of HH 158 and knot E of HH 702 in PVDs : [S n] A6731 
/ A6716, [N ii] A6548+A6584 / [O i] A6300+A6363, and 
[O i] A6300+A6363 / [S ii] A6716+A6731. The ratio 
of the [S Ii] doubl et can be use d to estimate the elec¬ 
tron density (lOsterbrockl Il989th which is larger than 
2 within ±5" and decreases with the distance. The 
[N ii] /[O i] ratio, which is proportional to the ion¬ 
ization fraction of hydrogen, is less than 0.5 at the 
center and increases at distances further than ~ 5". 
The [O l]/[S ii] ratio shows a similar trend with the 



Figure 8. Line ratios estimated from HH 158 and HH 702, 
for four distances: d = 0”, 7”, 14", 650". 

[S ii] doublet, large near the center and less than 1.0 
further than 5". This ratio is related to the electron 
temperature. The line ratios of HH 702 are compara¬ 
ble with those obtained from HH 158 at ~ 14" from 
the source. In Figure [S] we show the relations between 
three ratios, [S ii] A6731 / [O i] A6300, [N ii] A6584 
/ [O i] A6300, and [S ii] A6716 / A6731. In the fig- 
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Table 2 

Locations of the knots with respect to the source position during the last 30 years. 


Observation date 

Bl a 

Distance (") 
B2 a 

C D 

Reference 

1984 Oct. 9 

- 

- 

8.5 

Mundt et al. (1987) c 

1986 Dec. 

- 

- 

8.7 

Eisloffel & Mundt (1998) 

1992 Oct. 

2.25 

3.25 

- 

Solf & Bohm (1993) 

1994 Nov. 3 

2.7 

3.9 

- 

Lavallcv et al. (1997) 

1997 Jan. 

3.3 

- 

- 

Dougados et al. (2000) 

1998 Jan. 23-26 

3.6 

- 

- 

Lavallev-Foiinnet et al. (2000) 

1998 Feb. 5 

- 

- 

10.65 

HST d 

1999 Jan. 14 

3.8 

- 

- 

Maurri et al. (2014) 

2002 Sep. 17 

4.6 

- 

12 

Whelan et al. (2004) 

2012 Nov. 15 

6.74 

8.75 

13.7 

This work 

2014 Nov. 21 

7.65 

- 

13.5 

This work 

Proper motion ("yr _1 ) 

0.230 

0.272 

0.180 

- 

Ejection date 

1981 

1980 

1939 

- 


a knot name from ILavallev et all 1119971 1. 
b knot name from lEisloffel Hr. Mundtl 1199^1 . 

c position of knot C measured from Figure 9 of IMundt et all 119871 1. 

d data taken from the HST archive, image of WFPC2 with the F675WF filter (proposal ID. 6855). 


ure, the four points represent the values at distances 
d = 0", 7", 14", 650" from the source. The distance 
of 650" is the location of the knot E of HH 702. The 
ratios i ncrease with distance , as in a previous study 
(lLavallev-Fouauet et alJl2000l Figure 3a-b). 


4. Discussion 
4.1. Velocity Variation 


In Figure [3j the HVC is getting slower with distance. 
The velocity decreasing rate derived in Section [3] is ~ 7 
km/". If the jet is continuously slowing, then it should 
stop at 35" from DG Tau. However, the knots of HH 
702 located at 650" show a radial velocity of ~ —80 
km s” 1 . 


Radial velocity ( v rad ) measurements of knot C in 
HH 158 were reported in s evera l previous stu dies. The 
knot was observed at 8"5 (IMundt et al.lll987l) and 8"7 
( Eisloffel fc Mundtlll998l) from the source in 1984 and 
1987, respectively. We converted the heliocentric veloc- 
iti es to LSR velo cities f or the comparison with our data. 
In iMundt, et al.l rfl987h . Vrad were —136 k m s 1 in Ha 
and —123 km s -1 in [S n] A6731. In lEisloffel fe Mundt 
(119981) . [S ii] is ~ —146 km s” 1 . In 2002 ( Whelan et al. 
120041) . [S n] A6731 and [O i] A6300 show v rad of ~ -150 
and —200 km s -1 , respectively. For the three emission 
lines mentioned above (Ha, [S ii], [Oi]), our data show 
Vrad of —190, —128, and —180 km s -1 , respectively. 
Since different emission lines originate from different 
excitation conditions (i.e., temperature, density, etc.), 
they represent various regions in the outflow and show 
peaks in different velocities. The temporal variation of 
the velocity indicates that the velocity of knot is vari¬ 
able not only along the distance but also in time. This 
may imply that the ejected knot has been decelerated 
when it passes through the ambient medium. The ve¬ 
locity change of the knot could also occur due to the 
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Figure 9. Positional variations of three knots. For knot 
Bl, the data points are expanded from those in 'Pvo et aD 
(1200311 . 


shock caused by the interaction between the faster gas 
stream behind and the slower front knot. 


4.2. Proper Motion 

In the emission lines of HH 158, we have distinguished 
three knots at 6" 74 and 8" 75, and 13'/7 from the source 
in our data taken on Novem ber 1 5, 2 012. We identified 
these kno ts as Bl and B2 oflLavallev et ahl (|l997), and 
knot C of lEisloffel fe Mundtl ( 19981) . Recently. we could 
measure the knot positions from our additional data 
taken on November 21, 2014. For knot Bl and C, it 
was measured to 7'/65 and 13'/5, respectively. We note 
that knot B2 is not resolved in our spectra in 2014. 
The distance of knot C from the source is 0'/2 smaller 
in 2014 than that in 2012. The seeing size is > 2", so 
it may be due to the uncertainty of the measurement. 
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Figure 10. Flow chart showing how we derive the electron 
density (n e ), ionization fraction (a; e ), and electron temper¬ 
ature ( T e ). 


In Table [2] we list the locations of the knots dur¬ 
ing the last ~ 30 years including our work. Basi¬ 
cally, the posit i ons of knot B1 are extended from Ta¬ 
ble 1 of iPvo et aD ( 2003lh In the case of kn ot C, 
the f irst detection was reported in 1984 (iMundt et al.l 
1 19871) and it wa s seen at 12" in long-slit spectroscopy 
of AVhelan et al.l ( 2004 1 whose data were taken in 2002. 
In addition to these data points, we found an Hubble 
Space Telescope (HST) archival image obtained in 1998. 
The central peak in intensity contour of the knot in the 
HST image is located about 10"65 from DG Tau. 

Figure [9] shows the temporal variation of the distance 
for the knots Bl, B2, and C. We estimated the proper 
motion by linear regression: 0"230 ± 0"006 yr” 1 for 
Bl, 0"272 ± 0"002 yr' 1 for B2, 0"180 ± 0"008 yr" 1 for 
C (see Table [2]) . For knot Bl, the estimated value is 
smaller than t hose o btained from previous studies, e.g., 
0"272 in iPvo et akl (120031 ). Due to the large spatial 
resolution of our observation, the emission from knot 
Bl could be combined with that from knot Al which 
has identified from previou s studies dBacciotti et al.l 
1200(1 iDougados et al.ll2000l) . The smaller proper mo¬ 
tion value of knot Bl ma y be esti ma ted be cause knot 
Al is located 2"5 closer (jMaurri et al.1120141) to the star 
than knot Bl. From the radial and the tangential ve¬ 
locities of each knots, the inclination angles ( i ) of Bl, 
B2, and C are calculated as 32?5 ± 0?7, 37?0 ± 0?2, 
32?3 ± 1?1 wit h respect to the lin e of sight. 

For HH 702, IMcGroartv et all (120071) estimated the 
tangential velocities of knots A and E as 129 ± (10 
— 28) and 186 ± (10 — 28) km s” 1 in the Ha emis¬ 
sion, which correspond to 0"19 ± 0"03 and 0"28 ± 0"03 
yr -1 in proper motion, respectively. From the radial 
velocity of ~ —80 km s” 1 and the tangential veloci¬ 
ties, the i of knot A and E are estimated as ~58?2 ± 
1?1 and ~66?7 ± 1?1 with respect to the line of sight. 
The i of HH 702 are ~21—35° different from those of 
HH 158. We cannot avoid the possibility that the two 
HH objects were elected from different sources although 
IMcGroartv et al.l ( 20071) postulated that DG Tau seems 
the most probable driving source of HH 702 based on 
the proper motion direction. The different direction or 


the curved sha pe o f a jet may indicate the precession 
of the jet dRaga et al.ll200li) . 

4.3. Line Ratios and Physical Parameters 


iMaurri et al.1 (120141) shows the line ratios and phys¬ 
ical parameters within 5" of DG Tau. Their 
[S ii] A6731/A6716 ratio is simil ar to o ur value. Knots 
B0 and Bl at 3"3 and 3" 8 in iMaurri et al.l (120141 ) cor¬ 
respond to Bl at 6"7 in our data, after considering the 
proper motion. For this knot, the [S ii] doublet ratio in 
our data is slightly smaller than t heir estim ation. For 
the [N ii] / [O i] ratio, the value in IMaurri et al.l ( 20141) 
is 0.01 — 2.5, in a similar range as our data. 

In Figure flT)l the flow chart showing how we estimate 
the electron density (n e ), ionization fraction ( x e ), and 
electron temperature ( T e ) is displayed. Figure QT] shows 
the derived physical parameters for HH 158 and HH 
702. The estimated properties of each knot show good 
agreement with the planar shock mo del (lHartigan et al.l 
1 19941 : lLavallev-Fouauet et, al.ll2000l) . For the variation 
of the jet with the distance (i.e., variation with time), 
the model that treats the time variabi lity of the jet 
( Raga et~aklll990HRaga fc Kofmanll 19921) could explain 
our result with the internal shock caused by a variable 
(or periodic) jet ejection. 

In our data, it is hard to distinguish the two veloc¬ 
ity components within ± 5" from DG Tau due to the 
limitations in spatial and velocity resolution. The LVC 
is usually dominant close to the source, while HVC is 
extended further as shown in Figure [3j With this point 
of view, it is noticeable that the x e of HVC at |Y| > 
4" is higher than that of the LVC at |Y| < 4" but n e is 
lower in HVC in Figure |TT| One_the cither ha nd, prev i¬ 
ous s tudies (lLavallev-Fouauet et al.ll2000l : IMaurri et al.l 
l2014h show that the n e and x e are larger in their high 
and medi um velocity interval than those in low velocity. 
iHamannl (Il994l) also showed that the two velocity com¬ 
ponents of 20 TTS and 12 HAeBe stars are different in 
T e and n e . These results indicate that the two velocity 
components have different physical conditions. They 
are supposed to originate from different parts of star- 
disk systems or have different launching mecha nisms 
dKwan fe Tademarulll995l : IPvo et aDl2003l 120061) . 

4.3.1. Electron density (n e ) 

We calculated the electron density using the typical 
method of [S ii] A6716/A6731 doublet ratio ( Osterbrockl 
ll989l) . This ratio is known to be insensitive to tempera¬ 
ture variation, so we assumed the electron temperature 
as ~ 10,000 K to estimate the electron density (n e ). In 
the bottom panel of Figure [T2l n e is shown at four dif¬ 
ferent distance positions. At the proximity of the star, 
n e is ~ 10 4 cm' 3 and decreases to ~ 500 cm” 3 at 14" 
from the source. In HH 702 at 650", n e is estimated to 
be ~ 200 cm” 3 . 


4.3.2. Ionization fraction (x e ) 

We esti mated the ioniz a tion f raction from the model 
grid of lHartigan et al.l ( 1994 )). which shows the re- 
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Figure 11. (Left) Electron density (n e ), (center) ionization fraction (x e ), and (right) electron temperature (T e ) of HH 158 
and HH 702. 


x e , T e , n e , n H of DG Tau Jet 



Distance [arcsecond] 

Figure 12. (From top to bottom) Ionization fraction (x e ), 
electron temperature (T e ), total hydrogen density (uh), and 
electron density (n e ) of DG Tau jet as a function of the 
distance from the source. 


lations between the ionization fraction and the line 
ratio of [O i] A6300/Ha, [S n] A6716+A6731/Ha, 
[N n] A6548/[0 i] A6300 in Figures 3, 4, and 6. We 
assumed the magnetic strength of 100 /iG, close to 
the average value in the model. The middle panel of 
Figure |TT] and the top panel of [12] show the ioniza¬ 
tion fraction obtained from the [N ll]/[0 i] ratio. Fig¬ 
ure m shows the x e values derived from three differ¬ 
ent ratios of [O i] A6300/Ha, [S n] A6716+A6731/Ha, 
[N n] A6548/[0 i] A6300. The x e values deduced by 
[O i]/Ha and [S Il]/Ha ratios are in the range of 0.02 — 
0.2 which are much different from the values deduced 
from [N li]/[0 i]. Since the Ha emission flux comes not 
only f rom hydrogen recombination but also from col¬ 
lision dBacciotti fc Kisldffell 1 1 OOOt iHartisraiil l2003lh the 
use of the Ha emission in the estimation of the ion¬ 
ization fraction may cause uncertainty. Thus, we used 
only the [N n]/[O i] ratio in the calculation of x e . 


In Figures (TTj and [12] x e varies from ~ 10" 3 to ~ 
0.5 with distance from the sour ce. According to Fig¬ 
ure 16 of lHartigan et all (199.4) which shows the re¬ 
lation between ionization fraction and shock velocity, 
the shock velocity surpasses 100 km s _1 with x e of ~ 
0.5 aroun d knot C , which has ~190 km s _1 radial ve¬ 
locity. lLavallev-Fouauet et al.l (2000) showed that x e 
was increasing from 0.02 to 0.6 in HVC and from 0.004 
to O.J_ in LVC within 4" from the source. Recently, 
iMaurri et al.1 (2014) reported similar results within 5" 
using HST/STIS data obtained in 1999. 


For HH 702, x e is greater than 0.2, which is similar 
to the value at 14" from the source. 
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[O i] emission line with the following equation: 
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Figure 13. Ionization fraction (x e ) derived from the line ra¬ 
tios of [ O il/Hq, [S if]/Ha, an d [N Il]/[0 i] based on the 
model of lHartigan et all (119941 1. 


4.3.3. Electron temperature (T e ) 

Electron temperature (T e ) values are shown in Figures 
ED and m a nd are estimated from the model grid of 
Figure 3.6 in lAgra-Amboagel (2009) which shows the 
relation between emission line ratios and electron den¬ 
sity, temperature, and ionization fraction. Using the 
x e values obtained from 14.3.21 we have calculated the 
temperature from the model. T e is higher than 15,000 
K in the central region and decreases to ~ 5,000 K at 
14". At HH 702, the temperature was estimated to be 
slightly higher than that measured at knot C of HH 
158. 


4.4. Mass-loss Rate 

The mass-loss rate o f DG Tau was e stima ted in sev¬ 
eral previous studies. lHartigan et al.l (1992) estimated 
the value of 3.2 x 10~ 7 Mm yr -1 from the [O i) A6300 
emission line flux. In lLavallev-Fouauet et al.l ( 200 0). 
they have obtained the total value of 1.4 x 10~ 8 M 0 
yr^ 1 from the low, intermediate, and high velocity in¬ 
terval. The lowe r limit of mass-loss rate estimated by 
iBeck et al . ( 2010lf using the high -velocity Br y ernis- 
sion was 1.2 x 10 -8 Mq yr -1 . lAgra-Amboage et al.l 
( 201 11 has estimated a total mass-loss rate from —300 
km s -1 to —50 km s -1 comp onents of the [Fe ill emis¬ 
sion to 3.3 x 10 -8 Mq yr -1 . iMaurri et al.l (12014H also 
has estimated a mass-loss rate in the range of 10 -8 to 
10 -9 Mq yr -1 within the region 0"7 from the source. 

To estim ate the mass-loss rate , we have adopted the 
method of lHartigan et al.1 ( 1995 b We can calculate the 
mass-loss rate with the equation below: 


Mj = MtotV±/Ia 


(1) 


V± is the projected velocity on the plane of the sky and 
Z_l is the projected jet size of the aperture. The total 
mass of the jet M^ot is estimated from the forbidden 


M T ot = 9.61 x 10 


-6 


i + 2s 

Up. 


L 6300 

L e 


Mq (2) 


We estimate the mass-loss rate close to the star in 
the region within ± 5'/4 (^750 AU) from the star. 
The equivalent width of [O i] A6300 emission line is 
measured to be 21 A. We assume R = 8.74 mag¬ 
nitude, which_is_thedereddened R magnitude of DG 
Tau in lHartigan et all (Il995l) . The line luminosity of 
the [O i] A6300 (L 630 o) converted into solar units is 
9.2 x 10- 3 L 0 . V± is estimated to be 75 km s 1 
for the measured radial velocity 120 km s -1 , apply¬ 
ing the inclination angle of 32° with respect to the 
line of sight. For l±, we adopted 6.07 x 10 10 km as 
the projected jet aperture size which corresponds to 
our slit-width of 2"9 at a 140 pc distance to the Tau¬ 
rus molecular cloud. The F WHM of the jet measured 
from iDougados et al.l (|2000|) is ~0 / . / 21—1"25 at a dis¬ 
tance of ~0 , /4—4" from the star. lAgra-Amboage et al.l 
(1201111 obtained the value of Ml" 1—0"8 within ~1"4 
from the source, for the medium velocity component of 
the blueshifted jet. The jet width increases with dis¬ 
tance in the vicinity of the source. We assume that the 
jet width in our estimation would be larger than F'25 
because we estimate a mass-loss rate within ± 5"4. l± 
of 2"9 should be regarded as an upper limit. The tem¬ 
perature is assumed to be 82 00 K, the average tem per- 
ature for the [O i] A6300 line ( Hartigan et al.)ll994h and 
the critical density n c in this temperature is 1.97 x 10 6 


cm 3 . For the electron density n e , 7 x 10 4 '’ TTl 3 


cm 


is 


used. The mass-loss rate, Mj , for the [O i] emission 
line is ~ 1.0 x 10 -7 Mq yr -1 . 

The Mj can be also estimated using the [S n] A6731 
emission line. Mtot is given by: 


M T ot = 1-43 x 10 3 ^ ^ ^ M 0 


(3) 


The equivalent width of the [S Ii] A6731 line is ~ 
4.5 A. The Mj from [S ii] is 1.1 x 10 -7 Mq yr -1 , 
which is almost the same as that of [O i] A6300. 
The mass-loss rates of [O i] A6300 and [S ii] A6731 
lines are a factor of 3 lower than the value of 
lHartigan et al.l ( 1995ll and is 3 — 7 times l arger than 
those estimated by Lavalley-Fouauet et al.l (|2000U and 


lAgra-Amboage et al. (2011). The calculated rate can 


differ with time because the measured line intensi¬ 
ties and velocities of the jet could change due to the 
variability of the accretion/outflow rate of YSO. The 
rate becomes about two times larger if we replace the 
jet aper t ure siz e in our calculation with the value in 
lHartigan et al.l (ll995lh which is 1"5. The difference in 
spatial reso lution between ours a nd higher-resolution 
(e.g., lAgra-Amboage et al.l l201l!) observations cloud 
cause the difference in the mass-loss rate. 

5. Summary 

We have investigated the kinematics and physical char¬ 
acteristics of HH 158 and HH 702 from the optical long 
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slit spectroscopy. From the analysis through the PVDs 
and line profiles depending on the distance from the 
source, the peak velocity varies in the range of —50 to 
—250 km s -1 and the line width was measured to 150 
to 250 km s _1 . The proper motions were measured for 
knot Bl, B2 and C of HH 158: 0"230, 0"272, and 0"180 
yr -1 for the knots, respectively. 

From our analysis of the line ratios, we find that the 
physical properties (n e , x e , T e ) are changing along the 
jet from the proximity to the star to knot C at 14": 
(1) n e varies from an order of 10 4 to 10 2 cm' 3 , (2) x e 
varies from ~ 0 to 0.4 at 14", (3) T e varies from 15,000 
K to ~ 5,000 K. It is noticeable that the line ratios 
and physical parameters of HH 702 do not show much 
difference from those obtained at knot C of HH 158. 
The mass-loss rate (Mj) is estimated to be ~ 10~ 7 M 0 
yr , similar to values obtained in previous studies. 
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